Vibrational effects in laser driven molecular wires 
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The influence of an electron-vibrational coupling on the laser control of electron transport through 
a molecular wire that is attached to several electronic leads is investigated. These molecular vibra- 
tional modes induce an effective electron-electron interaction. In the regime where the wire electrons 
couple weakly to both the external leads and the vibrational modes, we derive within a Hartree-Fock 
approximation a nonlinear set of quantum kinetic equations. The quantum kinetic theory is then 
used to evaluate the laser driven, time-averaged electron current through the wire-leads contacts. 
This novel formalism is applied to two archetypical situations in the presence of electron-vibrational 
effects, namely, (i) the generation of a ratchet or pump current in a symmetrical molecule by a 
harmonic mixing field and (ii) the laser switching of the current through the molecule. 

PACS numbers: 05.60.Gg, 85.65.-|-h, 05.40.-a, 72.40.-l-w 



I. INTRODUCTION 

In recent years, considerable experimental progress in 
the determination the current-voltage characteristics of 
molecular wires has been achievediS*^ In these exper- 
iments, a single molecule is contacted with two nano- 
electrodes such that a transport voltage can be applied. 
One recent measurement, for example, focused on the 
influence of the chemical anchor group which couples 
the molecule, an oligothiophene derivative, to the elec- 
trodei^ In another experiment, the length-dependence of 
the current through a DNA strand^ was investigated. 
Most descriptions of such transport experiments rely on 
generalizations of the scattering approach put forward 
by Landaueri^ Presently, the main theoretical focus 
lies on the ab-initio computation of the orbitals relevant 
for the motion of excess charges through the molecular 
wire i2ii5:iM?i^ 

Another line of research employs rather generic 
models to gain a qualitative understanding of the 
transport mechanisms involved. An important prob- 
lem addressed in this way is the conduction mech- 
anism in the presence of electron-vibrational cou- 
p^^g,i4,i5,i6,i7,i8,i9,2Q,2i,22,23,24 ^j^j^ increasing strength 

of the coupling between the wire electrons and the vi- 
brations, the electrons tend to localize on single wire 
units. Correspondingly the transport mechanism changes 
from a purely coherent transfer to a sequential hopping 
process. The related transfer rates are known from the 
theory of non-adiabatic electron transfernS^ Such a spe- 
cial regime of charge transmission can also be described 
within a Redfield theory after carrying out a so-called 
polaron transformationi22i2& If the wire has to be de- 
scribed by spatially extended molecular orbitals, sequen- 
tial transfer proceeds by jumps from one electrode into 
a wire orbital, possibly followed by an intra-wire re- 
laxation, and then by a jump from a wire orbital into 
the other electrode. In contrast, the so-called superex- 
change mechanism enables a direct transfer from one elec- 
trode to the other. The crossover from superexchange 



to a sequential transfer mechanism has been studied in 
Refs. ll4l2ll23l24l27l2Sl In that context, a treatment 
of inelastic scattering processes profits from an unified 
theory of the electron transfer through molecular donor- 
acceptor complexes i2S 

Molecular wires illuminated by laser fields have been 
proposed for the investigation of ac transport effects like 
coherent quantum ratchetsj^S*^ optical control of current 
and noise;22i2^ and resonant current amplificationi^^iS^ 
The appropriate treatment of these ac phenomena is 
based on Floquet theory, which allows to take into ac- 
count the action of the time-dependent field exactly. A 
Floquet scattering approach for the fully coherent trans- 
port regime has been developed, '^'^ but it cannot be gen- 
eralized straightforwardly to the case with additional 
electron-vibrational coupling. Better suited for this sit- 
uation is a quantum kinetic equation formalism which, 
however, is perturbative in both the wire-lead coupling 
and the electron-vibrational coupling. 

In the absence of an external transport voltage, a driv- 
ing field can induce a so-called pump or ratchet cur- 
j.gj-^|.|36j37j38j39 rpj^g same happens even in perfectly sym- 
metric conductors if one adds a higher harmonic to the 
driving field. The investigation of the corresponding ef- 
fect for the motion of a particle in a tight-binding lattice 
revealed that the resulting current depends sensitively 
on the relative phase between the two components of the 
driving. For this system, two limits have been studied: 
the fully coherent dynamics and the overdamped Brow- 
nian motion. The dependence of the current on the rela- 
tive phase is in each case qualitatively different i^^iii The 
present model has the advantage that it enables the study 
of such effects also in the crossover regime between the 
purely coherent and the purely incoherent transport. 

A further intriguing phenomenon in driven transport is 
the suppression of the dc current caused by properly tai- 
lored ac fields. This effect is the transport counterpart 
of the so-called coherent destruction of tunneling (CDT) 
found in bistable potentials without any connection to 
external Ieads42ii2il4 When coupling the bistable system 
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to a heat bath, tunneling becomes a transient, see Ref.'45! 
for a review. The same is true for coherent destruction 
of tunnehng: ultimately, the driving-induced localization 
decays via dissipative transitions.''^ Here, we address the 
role of dissipation for the corresponding transport effect. 

The paper is organized as follows: In Sect. ^ we in- 
troduce a model for the laser driven molecule coupled to 
several leads and vibrational degrees of freedom. Sub- 
sequently, in Sect. IIIII we derive a kinetic equation ap- 
proach in the Floquet basis. An expression for the result- 
ing current through the molecule is derived in Sect. IIVI 
Finally, in Sect. the formalism is applied to study 
the influence of the vibrational coupling on non-adiabatic 
pumping and coherent current suppression. 



II. THE MODEL 

In the following, we consider a molecular wire which is 
attached to a number of electrodes and which is driven by 
an externally applied ac field. We neglect hole transport 
and start with the Hamiltonian 

Ho{R) = Eo + Hvih , (1) 

of the neutral wire, i.e., in the absence of excess elec- 
trons injected via the electrodes. Here, R denotes the 
set of all involved vibrational coordinates with equilib- 
rium configuration Rq and Eq = iJo(-Ro) is the electronic 
ground-state energy of the neutral wire (set equal to zero 
in the following). In a representation by normal modes 
with mode index ^, we find the vibrational Hamiltonian 

H^ih = ^ hLu^ blb^ (2) 

with the usual harmonic oscillator operators bl and b^. 
The normal modes may be delocalized over the whole 
wire or may be localized on specific wire units. We as- 
sume that they always remain in thermal equilibrium and 
are, thus, characterized by the Bose distribution tib (hujf^ ) . 

We describe the presence of an excess electron on the 
wire in the representation of N localized wire orbitals 
n), n = 1,...,N, and the corresponding Hamiltonian 
Hnn'{R,t)- In the absence of the external driving, the 
eigenstates of -ff„„' are the so-called LUMOs (lowest un- 
occupied molecular orbitals). To account for electron- 
vibrational coupling this quantity is expanded with re- 
spect to deviations from the vibrational equilibrium con- 
figuration Rq. To lowest order this results in 

ih+bl) (3) 

The matrix Hnn' (t) = Hnn' {Ro , t) taken at the vibra- 
tional equilibrium configuration Rq includes the inter- 
action with the external field. The T-periodic time- 
dependence Hnn' (t) = Hnn' {t + T) modcls the action 



of the external field on the excess electrons when mov- 
ing through the wire. We have in mind a dipole-type 
coupling between electron and field, i.e., a contribution 
eE{t)xnSn,n' to the Hamiltonian Hnn'{t), where E(t) is 
the electric field strength, x„ the position of site n, and 
— e the electron charge. The considered action of an ex- 
ternal field may in principle also induce a heating of the 
electron gas in the metallic leads as well as a thermal 
expansion of the leads. Of course, this can be of im- 
portance for any experimental realization of the studied 
mechanism. Nevertheless, the present approach will not 
attempt to account for such effects. 

We assume that the whole set of vibrational coor- 
dinates discerns into subsets, labelled by an index z/. 
Therefore, we write the electron-vibrational coupling in 
Eq. ijSJl as 

J2 (k + bl)=J2 J2 ^^-i + ' (4) 

where the notation ^ G indicates that the summation 
runs only over the modes in the corresponding subset. 
The related spectral densities read 

-D.(a.) = J^lAf.^p ,5(^-^4). (5) 

For notational convenience, we define the spectral density 
for w < by D^{—uj) = —D^iuj). In our numerical 
calculations, we restrict ourselves to a situation where 
each orbital \n) couples to exactly one of these subsets, 
i.e., 

-^nn'iy — ^nn' ^nu : (6) 

and assume identical Ohmic spectral densities D^ito) = 
Khuj, where the dimensionless coupling strength k is iden- 
tical for all sites. This model has been employed recently 
for the description of dephasing and relaxation in (time- 
independent) bridged molecular wires 

In order to take the exclusion principle properly into 
account, we employ a many-electron description based 
on a second quantized notation. Neglecting the Coulomb 
interaction among the electrons, the wire electrons are 
described by the Hamiltonian 

n,n' 

The operators c„ and cj^ annihilate and create, respec- 
tively, an electron in the orbital \n), n = 1,. . . ,N. We 
presume vanishing overlap among different orbitals such 
that the annihilation and creation operators obey the 
standard anti-commutation relations [c„,c|j,]+ = 6nn'- 
The electron-vibrational coupling Hamiltonian assumes 
the form 

Hel-vih = ^""'^ + 4) ■ 
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Next, we consider the coupling of the wire to L differ- 
ent macroscopic electronic leads described by the Hamil- 
tonian 

H\eads — ^ ek£ Ck£ ■ (9) 

The electrons in each lead t are labelled by wavevectors 
k referring to bulk or surface states with energy ek^. All 
lead states are mutually orthogonal (and also orthogonal 
to the wire states) and, therefore, the creation and an- 
nihilation operators Ck£ and cj^^, respectively, obey the 
standard anti-commutation relations. We assume that 
the lead electrons stay in thermal equilibrium and are, 
thus, described by the Fermi distributions /(ek£ — M*") 
with a common temperature T but possibly different 
electro-chemical potentials /x^. 

The coupling of each lead to exactly one of the suitably 
labelled molecular orbitals is described by the Hamilto- 
nian 

-ffwirc-lcads = ^ Vkf 4f Q + h.C. (10) 

with Vk^ being the tunneling matrix elements. As it will 
turn out, the coupling to the leads is completely charac- 
terized by its spectral spectral density 

= f El^kd'<5(6-ek.) , (11) 

k 

which becomes a continuous function of the energy e if 
the lead states are dense. Since we are mainly interested 
in the behavior of the molecular wire itself and not in 
the details of the lead-wire coupling;^ we will assume 
for all numerical calculations that the conduction band 
width of the leads is much larger than all remaining rele- 
vant energy scales. In this so-called wide-band limit, the 
spectral densities are constants, r^(e) = Ti. 

The dynamics of the present model is now fully speci- 
fied by the Hamiltonian 

H{t) = i7wirc(^) + -H^vib + ^Icads + -f^ol-vib + -f^wirc-loads ■ 

(12) 

In the following, we start from the uncoupled subsystems 
^wiro(^)+^vib+^ioads and treat the influence of both the 
wire-lead coupling and the electron- vibrational coupling, 
subsumed in the Hamiltonian 

HcovL\y\ -^c\-vih ^" -^wire-leads i (-^'^) 

within a master equation approach in second order per- 
turbation theory. 

III. QUANTUM KINETIC EQUATION 
APPROACH 

Within this work, we focus on the case of weak and in- 
termediately strong coupling of the electrons on the wire 



to both the vibrational environment and the electronic 
states in the leads. Thus, it is favorable to choose a de- 
scription in terms of the reduced density operator g^i^c (t) 
of the wire electrons which follows from the total density 
operator by tracing out those degrees of freedom which 
correspond to vibrations and lead electrons. Then, the 
derivation of a closed equation for Pwire (t) to second order 
in -ffcoupi represents a standard procedure of dissipative 
quantum dynamics (for a review, see, e.g., Ref. 25). We 
neglect initial correlations between the wire electrons and 
the environmental degrees of freedom, which stay in ther- 
mal equilibrium, and do not account for non-Markovian 
dissipative effects. The resulting quantum master equa- 
tion thus reads 

f?wire (^) r[-f^wire(Oj £?wire(0] 

n 

OG 

Jdr Tronv[-ffeoupl, [Hcoupl{t-T, t), Qwirc{t)<^gcnv,cq]]- 


(14) 

The trace refers to all environmental states, i.e., the elec- 
tronic states of the leads as well as the wire vibrations, 
and the operator i?coupi(i, t') = UQ{t,t')HcoupiUo{t,t') 
describes the coupling Hamiltonian H13|) in the interac- 
tion representation. The related zeroth-order time evo- 
lution operator 

= Texp (^-i ^ d<" [i7wire(i") + i?vib + ifleads]) (15) 

is responsible for the dynamics of the uncoupled subsys- 
tems. In this way the external driving field not only de- 
termines the coherent dynamics of the wire electrons but 
also the dissipative part of the master equation (|14|l 

Note that Eq. (|14|l still determines the dynamics of the 
full many-particle density matrix of the wire electrons. 
Later on, we will derive an equation of motion for the 
single wire-electron density matrix defined as 

Pnn'{t) = Trcl [gwiic{t) cl, Cn] , (16) 

where the trace runs over the many-particle states of the 
wire electrons. It will be demonstrated below that in 
the presence of an electron- vibrational coupling, a closed 
equation for P„„' (t) can only be obtained when an ap- 
proximation is carried out for the two-electron density 
matrices by employing a decoupling scheme. We will 
detail this point when introducing the Hartree-Fock ap- 
proximation in Section IlII CI However, before doing so, 
we will briefly review the Floquet method for the treat- 
ment of the explicit time-dependence appearing in the 
propagator ifT^ . 
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A. Introduction of Floquet states 

For the evaluation of Eq. (|14|1 it is essential to use an 
exact expression for the zeroth-order time evolution op- 
erator Uo{t,t'). The use of any approximation bears the 
danger of generating artifacts, which, for instance, may 
lead to a violation of fundamental equilibrium proper- 
tiesiSSii2 In the present case, the only non-trivial contri- 
bution to the propagator l(T^ stems from the periodically 
time-dependent wire Hamiltonian (t) = _ff„„' (<-|-T). 
A proper tool for the efficient computation of the corre- 
sponding propagator is Floquet theory (^SiSLSiiS which 
is based on the fact that there exists a complete set of 
solutions of the form 

= e-'^"*/'' |<i>„(t)) , = + T)) 

(17) 

with the so-called quasienergies and corresponding 
Floquet modes The Floquet modes fulfill the 

eigenvalue equation 

(^|n}iJ„„,(0(n'l-i^^^)l*aW) = ea|$a(0) • (18) 

The practical usefulness of the Floquet ansatz H17|l is 
rooted in the fact that the Floquet modes are periodic 
functions of time t and can therefore be decomposed in 
a Fourier series 



l*a(0) =E 1*"^'')' 
k 



(19) 



<i>a(i)) 



This representation makes explicit that each quasi- 
energy is equivalent to the quasienergies 



khn. 



(20) 



where k is an arbitrary integer. Thus, we can restrict 
ourselves to states with eigenvalues in one Brillouin zone, 
E - hn/2 <ta<E + Kl/2. 

It is now convenient to define a "Floquet picture" via 
the time-dependent transformation of the fermionic cre- 
ation and annihilation operators 



The inverse transformation 



n|$a(t))c„(t) 



(21) 



(22) 



This is readily verified by differentiating expression H23() 
with respect to t and using the fact that \^a{t)) is a so- 
lution of the eigenvalue equation (|18|l . The proof is com- 
pleted by noting that Eq. (|23|l fulfills the initial condition 

£a(i',i') =Ca(i')- 

It is advantageous to change to the "Floquet repre- 
sentation" of the single (wire) electron density operator 



PcP (t) =J2^n I (t)) (t) I n') Pnn' (t) 

n,7l' 

= Trcl [£'wiro(i) 4(t) Ca{t)] = {cjj{t)Ca{t)) t . 



(24) 



After some algebra, we obtain for the dynamics of these 
expectation values the expression 



^Pa/3(t)= -^(ea-e^)Pa/3(t) 

oo 

~ ^ ydr([[c]3(i)cQ(i),i/coupi],^coupi(i-T,t)])t 



(25) 



Obviously, the canonical transformation H21I) to the basis 
of the Floquet operators Ca{t) has diagonalized the co- 
herent part of the master equation (|25|l and the only task 
left is the evaluation of the incoherent contribution. Here, 
we use the fact that the contributions resulting from the 
coupling of the wire electrons to the electronic leads and 
to the wire vibrations can be treated separately. This is 
possible due to the assumption that the lead electrons 
and the vibrations remain uncorrelated at all times. Fi- 
nally, we will obtain as a main result a quantum kinetic 
equation of the form 



Pa,3 - -T(ea-e/3)^'a/3+^'a/3|^i,,_i^^ds + ^"/3|cl-vib- (^6) 



Specific expressions for the second and third term on the 
right-hand side of the last equation will be derived in the 
following sections. 



follows from the mutual orthogonality and the complete- 
ness of the Floquet states at equal times. '^'^ Note that the 
right-hand side of Eq. (|22|l becomes i-independent after 
the summation over a. In the interaction picture, the 
operators Ca{t) obey 

c„(t,i') = uUt,t') Co.it') Uo{t,t') = e-'^°('-*')/''c„(t') . 

(23) 



B. Floquet-representation of the wire-leads 
coupling 



For the evaluation of the contribution of the wire-leads 
coupling to the kinetic equation (|26|) , we have to evaluate 
the integral in Eq. H25|) for the corresponding term in the 
coupling Hamiltonian (|13|l . Using the relation Cqi{t — 
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T,t) = exp(ieg^r/?i) Cqi, we obtain 



C. Floquet-representation of the electron 
vibrational coupling 



-Pa/3 



wirc-lcads 



;{c-'-/'-'[([[4(0c„W,4],Q] + ),/(e-M£) 



(27) 



r/h 



[[[cpit] Ca{t),ce],cl] + )t f{e - fii) 
-{4[clit) Co.it), ce])t]} . 



For notational compactness, the time arguments of the 
interaction picture operators ci{t — T,t) have been sup- 
pressed. 

Using the transformation (|22|l . the commutators in 
Eq. |[77|) are readily evaluated to read 

[cl{t)ca.{t),ci] = -(£|$/3(i)}c„(i) , (28) 
[4(0c„(0,4] = ($„(0K)4(i) . (29) 

Moreover, the transformation (|22|l yields in conjunction 
with Eqs. H19() . H20() . and (|23|l the spectral decomposition 
of the wire operators in the interaction picture, i.e.. 



Ci 



-ikQt Ji.e.a^kT /h 



ma,k) Cc.it). (30) 



ak 



With the aid of the last two equations, one may readily 
carry out the time and the energy integrations in Eq. H27|) 
to obtain 



Pa/3 



wirc-lcads 



■k)Qt 



= 1 kk' 



(31) 



x{re{ec.,k){<^>a.k'\i}{e\^0,k) fi 

+ T,iep,k){^c^k' \e){i\^i3,k} fiep,k - lie) 

-^r,(ea',fc)($a,fc'|^)(^|$a',fc)Pa'/3(i) 
a' 

- J2 ^fi^P',k'){<i>p',k'\e)mp.k) Pcp'it)} 



Here, principal value terms stemming from an energy 
renormalization due to the coupling to the leads have 
been neglected. The terms containing Fermi functions 
describe resonant tunneling of electrons from the leads 
onto the wire, while the reverse processes are captured 
by the terms proportional to Papit). 



For ease of notation, we introduce the "position" and 
the "force" operators 



n,n' 

F,^Y.^'hAK+bl) 



(32) 
(33) 



respectively. Then, the vibrational contribution to the 
dissipative kernel in Eq. H25|l becomes 

-^"/^Icl-vib ^ 

OQ 

^-(^) Cc.it),X,lX,{t - T, t)])t 



OQ 

'-Y.j dr A,(t) {[[clit) cc.it), X,lX,it - T,t)] 



(34) 



The symmetrized and anti-symmetrized autocorrelation 
functions 



S.ir) 



1 

2h 



([F,(r),F.] + )c 



(35) 



D^iu) coth(?i(jj/2fcBT') cos(ti;r) , 



^.W = ^([F,(r),F,])eq 



(36) 



/ — Oi,iLo) sm(a;T) , 



of the "force" operators fully characterize the fluctuation 
properties of the wire vibrations. Note that coth(a;) = 
1/x + Oix) such that for an Ohmic spectral density no 
infrared divergence occurs in the integral H35(l . 

For the further evaluation of Eq. H34|l , we express the 
operator X^ and its interaction picture version X^it — 
r, t) in terms of the Floquet picture operators Cq, it) at 
time obtaining 

E E ^a/3.fe 4(0cMt) , (37) 

aJ3 k 

X^t -r,t)=J2J2 gi(E3-^=-fc?iO)r/n 

a,p k 

xci(<)c/3(i). 



^al3,k (38) 



The time-averaged coupling matrix elements in the Flo- 
quet basis have been abbreviated as 

- E E<'^".'^-+^' \^)Xnn'An'\^p,k') ■ (39) 

n,n' k' 

When evaluating Eq. (|34|l . it turns out that in addi- 
tion to terms containing the single-electron density ma- 
trix Pafsit), two-electron expectation values of the form 
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{cl{t) cl^{t) c/3{t) Ca{t))t appear. By iteration, one thus 
generates a hierarchy of equations up to iV-electron ex- 
pectation values. To obtain a description in terms of 
only the single-electron expectation values, we employ 
the Hartree-Fock decoupling scheme defined by the ap- 
proximation 

(clit) clit) cp{t) c^{t})t « P^s{t}Ppj{t) - Pps{t)P^^{t) . 

(40) 

Clearly, such a mean-field approximation only covers cer- 
tain aspects of the full many-particle problem. Never- 
theless, it offers a description which is consistent with 
the second law of thermodynamics, as we will detail in 
Sect. IIII El We remark that in principle one could also 
include electron-electron interaction in the framework of 
the mean-field approximation H4U|) . similar to the ap- 



proach put forward in Refs. Is^ls^ . However, focusing 
on the vibration-mediated interaction effects, we here re- 
frain from doing so. 

For the description of the transport problem, a scat- 
tering approach, i.e., a strict one-particle picture, is fre- 
quently employedii^*i^*2ii2LSiS Then, non-linear terms 
of the type H40(l do not appear. The same happens if one 
considers a closed system with a single electron; then, 
the corresponding equation is also of the form f^I^ but 
without the terms quadratic in PQ,^(t)4LSLS 

Upon insertion of the Hartree-Fock approximation H4()|l 
into Eq. (|34|l and disregarding again the principal value 
integrals, which correspond to an energy renormalization 
due to the electron- vibrational coupling, we obtain after 
a straightforward calculation 



PcfS 



cl-vib 



a'P' 



a" 0" 

2 E [E^^""'^""*^^) ~ E {^*a"0"aa'i't) +^l3"a'a"a{t)^Pa"f3" 
a' 13" q"/3" 



Pa'l3' 



13' a" 



a"f3" 



Pal3' 



r 



(41) 



Here, we have introduced the time-dependent, complex- 
valued coefficients 



^a(3a'l3'{t) =2^51'^'''' "'^^^ ^af3 ,k' ^0' a' -k 



l,k' ^^0'a',-k 
V kk' (42) 

X N^iea -^13 + k'nn) , 



where the functions N^{e) are defined for each vibrational 
subsystem v by N^^e) = D^{e/h) n^{e)/h. For an Ohmic 
spectral density, Dy{e/h) oc e, Ni,{e) is well-defined in 
the limit e ^ 0. 



D. Rotating- wave approximation 

For a very weak coupling of the wire electrons to the 
environmental degrees of freedom, the coherent time- 
evolution dominates the dynamics (|26|l of the density 
matrix Paf3{t). More precisely, the largest time-scale 
of the coherent evolution, which is given by the small- 
est quasienergy difference, and the dissipative time-scale, 
which is of the order of the coupling rates r^(e), are well- 
separated, i.e.. 



(43) 



for all ^, fc, e and a ^ [3. Note that this condition is only 
satisfiable if the quasienergy spectrum has no degenera- 
cies. Then, it is possible to replace the T-periodic coef- 
ficients in Eq. (|26|l by their time-averages. Furthermore, 
we may assume that in the long-time limit, the solution 
becomes diagonal and time-independent, i.e., we make 
the ansatz 



Papit) = const. = Sap fa 



(44) 



With these approximations, the quantum kinetic equa- 
tion ^ assumes the form 



pleads 



+ E (1 - /") - E (1 - M • 

a' a' 

(45) 

Here, the population of the Floquet states a due to the 
wire-lead coupling is determined by the decay rates 



..leads 



= EEl(^l*"^^)l'r.(eo..) 



(46) 



«=1 k 
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and the source terms 

L 

pleads = ^ ^ I (£|ci,„ ,) |2 r,(e„,fe) - ^,). (47) 

1=1 k 

In addition, the elcctron-vibrational interaction con- 
tributes a quantum Bohzmann type coUision term to 
Eq. H45|) . which takes into account the Pauh principle 
by the blocking factors 1 — /q. The corresponding scat- 
tering rates from one state a' into another state a are 
given by 

= ^-(^" - + ^^") • (48) 

f k 



non-linear equations H26() . which however, is numerically 
rather time-consuming and, especially in the strongly 
driven regime, only applicable for not too large systems. 
We then verify the T-periodicity of the resulting solution 
and compute its Fourier decomposition 

Pa/3(<)-5I'^"'''*^«/3.fc- (52) 

As we will see in the subsequent section, the Fourier co- 
efficients 

Pa;3,fe-^y^ dte"="*P„;3(i) (53) 
fully specify the stationary current through the wire. 



E. Thermal equilibrium 



IV. CURRENT THROUGH THE WIRE 



An important consistency check of the present the- 
ory is a thermal equilibrium situation, where -ff„„' is 
time-independent and where no external bias is present 
{fie = /i for all £).^^ One can show that the diagonal 
ansatz (|44|l leads to an exact stationary solution of the 
kinetic equation l|26|) . The populations fa obey also the 
kinetic equation H45() . but with the Floquet states and 
quasienergies replaced by the (adiabatic wire) eigenstates 
\ipa) and the eigenenergies Ea, respectively, of the static 
Hamiltonian Moreover, only fc = contributes to 

the rates 

Thermal equilibrium with respect to the coupling to 
the vibrational subsystems is characterized by the de- 
tailed balance condition 

(1 - /a) /a' = <!l (1 - /a' ) /a , (49) 

where the rates w^'^, satisfy 



The net (incoming minus outgoing) current through 
contact £ is given by the negative time derivative of the 
electron number Ng = ^l^iCke in lead £ multiplied by 
the electron charge — e, 

I,it)^e^^{m),^e^{[Hit),m])^. (54) 

For the Hamiltonian (|12|l . the commutator in Eq. H54I) is 
readily evaluated to read 



= - Y 

k 



(55) 



For a weak wire-lead coupling, one can assume for all 
times a factorization of wire and lead degrees of freedom. 
This assumption allows one to derive from Eq. 1)55(1 an 
explicit expression for the stationary, time-dependent net 
electrical current through the contact £ in terms of one- 
particle expectation values of the wire electrons, 

oo 

Mt) = in^'' j j '^^ r,(e) e'"'^ [{c\ ce{t, t - t))^_^ 



-[4,ceit,t-T)] + f{e- lii)^ . 

(56) 

tional condition = s'°^<'^ that ^i' adjusts to the ^^te that the anti-commutator [4, c,(t, t- t)]+ is in fact 



,,vib 



^E^,-E^)/ksT 



(50) 



A solution of Eqs. ((49|l and 1)50(1 is given by the Fermi 
distribution fa — f{Ea — fJ-'), where the chemical po- 
tential /z' remains undetermined. It is due to the addi- 



chemical potential /i of the leads. Thus, the equilibrium 
solution of the kinetic equations l(2()(l reads 



fa = f{Ea - ^J.) , 



(51) 



in accordance with elementary principles of statistical 
physics of non-interacting fermions. 



F. Numerical solution of the kinetic equation 

For the solution of the nonlinear kinetic equation l(26() 
one generally has to resort to numerical methods. We 
therefore use a propagation scheme for the computa- 
tion of the long-time limit of the solutions of the set of 



a c-number, which by means of the transformation 1(22(1 
and the interaction picture dynamics l(23() of the wire 
operators in the Floquet picture reads 

[4,£,(t,t-r)]+=^c--°-/'"'($„(i-r)|^)(^l<i>„(0) . 

a 

(57) 

Similarly, the expectation value appearing in the current 
formula ((55(1 can be expressed in terms of the density- 
matrix elements ((24(1 as 

(4 c,{t,t- t))^_^ = ^ e---/'^ ($^(t - rmmait)) 

X Pal3{t - t) . 

(58) 



These relations together, together with the spectral de- 
compositions of the Floquet states and of the density 
matrix, Eqs. H19|l and (|52(l . respectively, allow one to 
carry out the time and energy integrals in expression (|56|l . 
The current Ii (t) obeys the time-periodicity of the driv- 
ing field. However, because we consider here excitations 
with frequencies in the optical or infrared spectral range, 
the only experimentally accessible quantity is the time- 
averaged current. In the wide-band limit — the extension 
to the general case is straightforward — we thus obtain 

ak I3k' (59) 
- |(^|*a.fc)|'/(ea,fc-/^£) . 

Charge conservation, of course, requires that the 
charge on the wire Qwiro(0 can only change by a current 
flow, amounting to the continuity equation Qwh-c{t) — 
J2e=i -^^i^)- Since the charge on the wire is bounded, 
the long-time average of Qwireit) must vanish. From the 
continuity equation one then finds 

L 

E^^ = 0- (60) 

For a two-terminal configuration, ^ = L, R, we can then 
introduce the time-averaged current 

/ = /l - -/r. (61) 

To close this section we consider the thermal equilib- 
rium situation described in Sect. IIII El It is character- 
ized by /if = /i (absence of any transport voltage) and 
\^a,k) = unless A: = (absence of driving). Inserting 
the equilibrium solution Paf3,k = Sap Sk,o f{Ea — m) into 
Eq. (|59|l , we immediately see that all currents Ii vanish — 
despite any possible asymmetry of the molecule itself or 
of its coupling to the environment. Thus, the results of 
our microscopic theory are in accordance with the second 
law of thermodynamics. 

V. CONTROL OF CURRENTS IN TWO-SITE 
SYSTEMS 

With the necessary formalism at hand, we study in this 
section two aspects of driven transport for which dissipa- 
tion plays a significant role: (i) the generation of pump 
currents by means of harmonic mixing fields and (ii) op- 
tical current switching. As a rather generic model, which 
still captures the essential physics, we employ a symmet- 
ric wire that consists of two sites in the two-terminal con- 
figuration sketched in Fig. ^ The wire sites are coupled 
by a hopping matrix element A while an electromagnetic 
field causes time-dependent shifts of the on-site energies. 
Then, the wire Hamiltonian reads 

-ffwire(t) = -A(C[CR + 4cl) + ^(c[cl " 4cr) , (62) 
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FIG. 1: Symmetric two-site structure coupled to two leads 
which rectifies an externally applied laser field of the form 
a{t) = Ai sm{flt) + A2 sm{2flt + if). 

where a{t) ^ a{t + T) represents the dipole force on the 
wire electron multiplied by the distance of the two wire 
sites. Furthermore, we assume that the molecule couples 
equally strong to both leads, thus, Fl = Fr = F. 

In a realistic wire molecule, A is of the order 0.1 eV. 
Thus, a wire-lead coupling strength F = 0.1 A/S corre- 
sponds to a current eF = 2.56 x 10"^ Ampere and a laser 
frequency = A/h lies in the infrared spectral range. 
Furthermore, for a distance of lA between two neighbor- 
ing sites, a driving amplitude A = A is equivalent to an 
electrical field strength of 10^ V/cm. 

It turns out that for the description of the effects dis- 
cussed below, the off-diagonal elements of the single- 
particle density matrix Pa/3{t), a ^ /3, play an essen- 
tial role. Therefore, we have to go beyond the rotating- 
wave ansatz (|44f) and consequently employ the propaga- 
tion scheme for the full nonlinear kinetic equation H26() 
described in Sect. IIII Fl 

A. Non-adiabatic pump current from harmonic 
mixing 

The Hamiltonian (j62|l with a driving of the form 
a{t) oc sin(f2t) has an intriguing symmetry, the so-called 
generalized parity, which has been discussed widely in 
the context of driven tunneling^^^iJi: a time translation 
by half a driving period results for the external field in 
a minus sign, i.e., a{t + T/2) — —a{t). Thus, for the 
dipole coupling given in the Hamiltonian (|62|l . the time 
shift t t + T/2 is equivalent to interchanging the left 
and the right wire site. In addition, the dc current / 
is also inverted. Consequently, because / has to be inde- 
pendent of any (finite) time translation, it must vanish.— 
However, for a time-dependent driving field of the form 

a{t) = Ai sin(m) + A2 sin(2f7i + ip), (63) 

with Ai,A2 7^ 0, the generalized parity is no longer 
present and a non-adiabatic pump current emerges from 
the harmonic mixing of the two driving frequencies^ 
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FIG. 2: Average current through the two-site wire from Fig.Q 
driven by the harmonic mixing signal 1)63^ with amplitudes 
j4i = 2A2 = A as a function of the wire-lead coupling strength 
r (for K = Q) for different values of the phase difference ip. The 
driving frequency is = A//i, and the temperature is k^T = 
0.25 A. The dotted line is proportional to F, corresponding 
to a current that is proportional to F^. 
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Its magnitude is generally proportional to the coupling 
strength between wire and lead, / oc F, with a prefactor 
that depends on the details like the phase lag ip or the 
amplitudes Ai, A2. The phase lag (p = represents a 
particular case because for this value the wire Hamilto- 
nian obeys time-reversal parity, i.e., it is invariant under 
the operation (L,R, t) — > (R, L,— t). As a consequence, 
one finds that the dc current vanishes to linear order in F 
such that / oc F^i^ Figure |21 demonstrates this behavior 
and, moreover, reveals that already small phase lags of 
the order cp « 0.001 are sufficient to alter the qualitative 
F-dependence of the dc current. 

Harmonic mixing has also been studied recently for 
the motion of a quantum particle in an infinitely ex- 
tended tight-binding lattice, both in the purely coher- 
ent regime^iiSS and for incoherent, sequential quantum 
between adjacent sites42i^ It turns out that the depen- 
dence of the current on the phase lag tp is qualitatively 
different in these two limiting cases. This raises the ques- 
tion how the phase dependence of the current changes as 
a function of the dissipation strength. 

Generally, quantum dissipation results from a coupling 
of the quantum system to an environment — here, the 
metallic leads and the vibrational modes. Figure 13^^ de- 
picts the influence of only the wire-lead coupling: For a 
very weak coupling strength F — O.OOIA//1 and k = we 
find a dc current proportional to sin p. With increasing 
F, the dependence on p shifts towards cos p. In order 
to investigate the influence of the vibrational coupling, 
we choose again F = 0.001 A/ft and a finite but small vi- 
brational coupling strength k. The (/3-dependence of the 
pump current is given in Fig. ^jp; it exhibits the same 
dissipation- induced shift. 

Interestingly enough, the electron- vibrational coupling 
can enhance the pumping effect. This enhancement is 
most pronounced in the presence of the time-reversal 
parity discussed above, i.e., for p = 0. Figure^ shows 



FIG. 3: Average current through the two-site wire sketched 
in Fig.0driven by the harmonic mixing signal 16311 as a func- 
tion of the phase difference p (a) for different wire-lead cou- 
pling strengths F (for k = 0) and (b) for different electron- 
vibrational coupling k (for F = O.OOlA/fi). All other param- 
eters are as in Fig. |5| 



2 


_ 1 1 1 1 

— F = 0.005A/fi 


1 


1 


— - F = O.OlA/?i 






^ F = 0.05A/?i / 




0.5 


; - - F = O.lA/ft 




0.2 






0.1 


" 1 1 1 1 


1 




lO""* 10"^ 10^2 


10-1 




K 





FIG. 4: Average current (in units of eF^) through the two- 
site wire sketched in Fig. Q driven by the harmonic mixing 
signal 16311 with amplitudes Ai = 2A2 = A and phase differ- 
ence (/3 = as a function of the electron-vibrational coupling 
strength k. Different values of the wire-lead coupling strength 
F are shown. 



the pump current as a function of the vibrational cou- 
pling strength k. We find that the dc current can be 
increased by more than one order of magnitude. For 
values K <C /iF/A, the main dissipation comes from the 
leads and the vibrations are practically without influence. 
Correspondingly, one is back to a the situation of Fig. |21 
where the pump current is proportional to F^. 
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FIG. 5: Two-site wire attached to leads with the chemical 
potential difference fin ~ fJ,h ~ eV. 



B. Laser-switched current gate 

An external driving field can not only induce a pump 
current through the molecular wire, but for proper pa- 
rameters can also cause the opposite effect: A driving 
of the shape a{t) = As'm{Ut) can suppress almost com- 
pletely the dc current even in the presence of a large 
transport voltage Vm^ The physics behind these sup- 
pressions is the so-called coherent destruction of tunnel- 
ing (CDT) that has been found in the context of tun- 
neling in time-dependent bistable potentials4242iiSiSiiSS 
The central phenomenon observed there is that for a driv- 
ing with amplitude and frequency such that the ratio 
A/hn equals a zero of the Bessel function Jq (i.e., for 
the values 2.405.., 5.520.., 8.654.., . . . ), the coherent tun- 
neling dynamics comes to a standstill. As a related ef- 
fect for the transport through such a tunnel system, one 
finds pronounced suppressions of the dc current^ The 
fact that coherent destruction of tunneling is disrupted 
by finite dissipation^^i^^iSiSiSM^ motivates our inves- 
tigation of the influence of dissipation on these current 
suppressions. 

We model the transport voltage V by shifting the 
chemical potential of the left (right) lead, /xl (a*r), by 
-eV/2 {+eV/2), cf. Fig. El Due to the external voltage, 
the electric field can in principle also cause a static bias 
to the wire levels. We do not take this effect into account 
in the present work, but remark only that in the absence 
of a vibrational coupling, the current suppressions are 
stable against an internal biasi^ 

Before focussing on the influence of electron- 
vibrational coupling, let us first substantiate the discus- 
sion of the current suppressions by the numerical results 
depicted in Fig.|^ The time-averaged current / as a func- 
tion of the laser amplitude A exhibits, besides a global 
decay, pronounced minima whenever the CDT condition 
is fulfilled, i.e., when the ratio A/Ml assumes a zero of 
the Bessel function Jq. However, the current does not 
vanish exactly, but a residual current remains; its value 
is proportional to the molecule-lead coupling F, as can be 




20 40 60 80 100 



FIG. 6: Average current (solid) vs. driving amplitude for the 
setup sketched in Fig. The leads' chemical potentials are 
/xr = —fiL = lOA; the other parameters read HQ = lOA, 
ksT = 0.25 A, fir = 0.1 A. The dashed line marks the result 
obtained within a rotating-wave approximation. The inset 
depicts the minimal current Jmin at the first suppression as a 
function of the wire-lead coupling strength F. 



appreciated from the inset of Fig. |S| Since the current in 
the undriven situation is also proportional to F, we thus 
can conclude that the maximal suppression is determined 
by a factor which is independent of F. 

The inspection of the quasienergy spectrum reveals 
that CDT is related to crossings of the quasienergiesi^ 
Thus, at the center of the current suppressions, the 
quasienergies are degenerate and the condition H43|l for 
the applicability of the rotating-wave approximation is 
not fulfilled. Indeed, the dashed line in Fig. demon- 
strates that a rotating-wave approximation fails com- 
pletely in the vicinity of current suppressions. 

A central question to be addressed is the robustness 
of the current suppressions against dissipation. In the 
corresponding tunneling problem, the CDT driving al- 
ters both the coherent and the dissipative time scale 
by the same factor..^ Thus, one might speculate that 
a vibrational coupling leaves the effect of the driving on 
the current qualitatively unchanged. Figure [3 however, 
demonstrates that this is not the case. With increas- 
ing dissipation strength k, the characteristic current sup- 
pressions become washed out until they finally disappear 
when K becomes of the order unity. This detracting influ- 
ence underlines the importance of quantum coherence for 
the observation of those current suppressions. Note that 
the dissipation affects only the depth of the suppressions 
while the width remains unchanged. We close this sec- 
tion with the remark that within the present setup of two 
driven tunnel- coupled orbitals (cf. Figure [SJ and within 
our chosen parameter range, we do not detect the ana- 
logue of the effect of a stabilization of CDT within a cer- 
tain temperature range^£iSS4SLSS4S2i2£ or, likewise, with 
increasing external noisej^i as it has been reported for 
driven, dissipative symmetric bistable systems. 
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A[A] 

FIG. 7: Average current vs. driving amplitude for the setup 
sketched in Fig. |5] in the presence of dissipation of the 
form @. The inset depicts the minimal current Jmin at the 
first suppression as a function of the electron-phonon coupling 
strength K for F = 0.001 A/fi (solid line), 0.01 A/fi (dashed), 
and 0.1 A/ft (dotted). All other parameters are as in Fig.|S| 

VI. CONCLUSIONS 

We have derived a nonlinear quantum kinetic equation 
that allows one to investigate for a molecule the simulta- 
neous influence of a laser field, a coupling to leads, and 
in addition, a coupling to vibrational modes. The use 
of Floquet states as a basis set for the reduced single- 
particle density matrix represents a most important tech- 
nical cornerstone. It enables one both the exact inclusion 
of the driving field and an efficient treatment of the dis- 
sipative couplings. Since the vibrational modes provide 



an effective electron- electron interaction, a formalism for 
general situations requires one to resort to further ap- 
proximations such as a Hartree-Fock decoupling scheme. 

Within this kinetic equation formalism, we have inves- 
tigated the influence of quantum dissipation on recently 
proposed transport effects caused by the action of laser 
flelds on molecular wires. For the non-adiabatic electron 
pumping that emerges from harmonic mixing, we flnd 
that dissipation can play a constructive role to the ex- 
tent that it can significantly enhance the current. 

For the model under investigation, we observed an en- 
hancement of the pump current by more than one order 
of magnitude. Moreover, the present scheme allows one 
to trace back the dependence of the pump current on the 
phase lag between the two harmonic mixing fields to the 
increasing influence of dissipation. The situation is less 
promising for effects that depend intrinsically on quan- 
tum coherence. We have found that the coherent current 
suppressions are derogated by the coupling to vibrational 
modes. Nevertheless, the effect persists provided that 
the quantum dynamics of the wire electrons remains pre- 
dominantly coherent. Finally, we share the hope that our 
general theoretical findings will provide motivation and 
prove useful to experimentalists to initiate correspond- 
ing, tailored experiments on driven molecular wires in 
the not too distant future. 
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